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Marine hydrocarbon seeps are ecosystems that are rich in methane, and, in some cases, 
short-chain (C2-C5) and longer alkanes. C2-C4 alkanes such as ethane, propane, and 
butane can be significant components of seeping fluids. Some sulfate-reducing microbes 
oxidize short-chain alkanes anaerobically, and may play an important role in both the 
competition for sulfate and the local carbon budget. To better understand the anaerobic 
oxidation of short-chain n-alkanes coupled with sulfate-reduction, hydrocarbon-rich 
sediments from the Gulf of Mexico (GoM) were amended with artificial, sulfate-replete 
seawater and one of four n-alkanes (C1-C4) then incubated under strict anaerobic 
conditions. Measured rates of alkane oxidation and sulfate reduction closely follow 
stoichiometric predictions that assume the complete oxidation of alkanes to CO2 (though 
other sinks for alkane carbon likely exist). Changes in the 8 13 C of all the alkanes in 
the reactors show enrichment over the course of the incubation, with the C3 and 
C4 incubations showing the greatest enrichment (4.4 and 4.5%o, respectively). The 
concurrent depletion in the 8 13 C of dissolved inorganic carbon (DIC) implies a transfer 
of carbon from the alkane to the DIC pool (-3.5 and —6.7% for C3 and C4 incubations, 
respectively). Microbial community analyses reveal that certain members of the class 
Deltaproteobacteria are selectively enriched as the incubations degrade C1-C4 alkanes. 
Phylogenetic analyses indicate that distinct phylotypes are enriched in the ethane reactors, 
while phylotypes in the propane and butane reactors align with previously identified C3-C4 
alkane-oxidizing sulfate-reducers. These data further constrain the potential influence of 
alkane oxidation on sulfate reduction rates (SRRs) in cold hydrocarbon-rich sediments, 
provide insight into their contribution to local carbon cycling, and illustrate the extent to 
which short-chain alkanes can serve as electron donors and govern microbial community 
composition and density. 

Keywords: short-chain alkanes, methane, ethane, propane, butane. Gulf of Mexico, microbial sulfate reduction, 
C1-C4 hydrocarbons 



INTRODUCTION 

Marine hydrocarbon seeps are natural features that support con- 
siderable biological diversity and activity (Muyzer and Van Der 
Kraan, 2008 and references therein). Tectonic activity forms faults 
that facilitate the release of methane (CH4) and other hydrocar- 
bons from deep subsurface oil and gas deposits, as well as gas 
hydrates (Muyzer and Van Der Kraan, 2008). Methane is a key 
constituent of the carbon cycle as it is one of the final products 
of the microbial degradation of organic matter (Thauer et al., 
2008 and references therein), though it can also be produced abi- 
otically through thermo chemical and geogenic reactions (Horita 
and Berndt, 1999). 

Because CH4 is a potent greenhouse gas, there is consid- 
erable interest in determining the fate of both biogenic and 
abiotic methane (for review see Reeburgh, 2007; Thauer et al., 
2008). Consequently, microbial methane oxidation under both 
aerobic and anaerobic conditions has received considerable atten- 
tion (Thauer et al., 2008 and references therein). The anaerobic 
oxidation of methane, or AOM, has been the subject of research 



for four decades, and much of the work has been focused on 
identifying the responsible microbes, their distribution, and the 
biochemistry of AOM. To date, five distinct mechanisms of AOM 
have been discovered (Callaghan, 2013 and references therein; 
Haroon et al., 2013). The AOM mechanism most relevant to 
this study is mediated by microbial consortia of archaea, related 
to the archaeal anaerobic methane oxidizer group ANME, and 
bacterial sulfate-reducers of the Desulfosarcinales/DesuZ/ococcus 
(DSS) group (Knittel and Boetius, 2009 and references therein). 
This microbial consortium mediates coupled AOM and sulfate 
reduction, though the exact nature of the association is not fully 
understood (Knittel and Boetius, 2009; Milucka et al., 2012 and 
references therein). 

The study of C2-C5 hydrocarbon degradation has lagged 
behind that of CH4 in spite of their abundance in some envi- 
ronments. Analogous to methane, C2-C5 gases are formed due 
to thermal cracking of fossilized organic deposits, and C1-C2 
gases are also produced biologically (Lorant and Behar, 2002; 
Hinrichs et al, 2006; Jones et al, 2008; Xie et al, 2013). At some 
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sites, including the Gulf of Cadiz, the Gulf of Mexico (GoM), the 
Caspian Sea, the Monterey Bay canyon (Lorenson et al., 2002), 
and the Norwegian continental shelf (Hovland and Thomsen, 
1997), C2-C5 hydrocarbons can account for 14-38% of the total 
gas (see Milkov, 2005 and references therein). In these areas, the 
oxidation of C2-C5 hydrocarbons may be a significant contrib- 
utor to the community bioenergetics of marine seeps (Lorenson 
et al, 2002; Formolo et al., 2004; Sassen et al., 2004; Alain et al, 
2006). For example, at GoM cold seeps, sulfate reduction rates 
(SRRs) are higher than can be accounted for by AOM alone, 
indicating that sulfate reduction is linked to the oxidation of 
other organic compounds potentially including short-chain alka- 
nes (Joye et al., 2004; Orcutt et al, 2010; Bowles et al, 2011). The 
extent of this as well as the influence of short-chain hydrocarbon 
oxidation on AOM is poorly constrained, but it is possible that 
C2-C5 hydrocarbon degradation is a significant process that co- 
occurs with AOM, and may compete for a common oxidant (i.e., 
SO^ ) (Joye et al, 2004; Orcutt et al, 2010; Bowles et al, 2011). 

Recent studies of marine and terrestrial seeps, as well as marine 
hydrothermal vents, have observed the microbial oxidation of 
short-chain alkanes coupled to sulfate reduction across a range 
of temperatures. The microbes, as revealed by phylogenetic anal- 
yses of isolates as well as enrichments, align with sulfate reducers 
within the Deltaproteobacteria and the Firmicutes (Kniemeyer 
et al, 2007; Savage et al, 2010; Adams et al, 2013; Jaekel et al., 
2013). Kniemeyer et al. isolated a bacterium, BuS5, allied to 
the DSS group within the Deltaproteobacteria that can anaer- 
obically oxidize propane and n-butane while reducing sulfate 
(Kniemeyer et al, 2007). Savage et al. showed that propane and 
n-pentane degrading enrichments from a terrestrial hydrocarbon 
seep were also dominated by the DSS group (Savage et al., 2010). 
Jaekel et al. further characterized propane and butane degrad- 
ing sediment-free enrichments to expand the understanding of 
the physiology of these microbes (Jaekel et al., 2013). Unlike 
these previous studies, Adams et al. observed appreciable rates of 
ethane degradation coupled to sulfate reduction in ex situ sedi- 
ment slurry-based batch reactors with sediments from the Middle 
Valley hydrothermal vent field (Adams et al., 2013). Collectively, 
these studies have placed some constraints on the relationship 
between alkane oxidation and sulfate reduction, though little 
remains known about the stoichiometric relationship between 
alkane oxidation — including methane- and sulfate reduction by 
mixed, natural communities, and their impact on local carbon 
cycling. 

The northern slope of the GoM is an ideal site to study the 
anaerobic consumption of short-chain alkanes because the sed- 
iments lie over hydrocarbon deposits including structure II and 
H gas hydrates rich in C1-C5 gases (Joye et al, 2004). The sites 
of hydrocarbon seepage in the GoM are also characterized by 
the presence of mats dominated by the chemoautotrophic sul- 
fur oxidizing bacterial genus Beggiatoa (Joye et al., 2004). The 
presence of these microbes suggests that H2S is available in the 
environment, which may indicate high advection and seepage 
rates (Joye et al., 2004). At the northern slope of the GoM, 
methane is the dominant component of seeping gas (72-96%) 
with some contribution from short-chain alkanes. The abun- 
dance of the short-chain alkanes decreases with chain length 



with ethane comprising 2.4-12.4% of the total gas, followed by 
propane and butane (iC4 + nC4) (1.2-12.6 and 0.3-4.3%, respec- 
tively) (Sassen et al., 1998). Stable carbon isotopic properties of 
the starting materials (vent gas from the deep subsurface), inter- 
mediate products (in situ gas hydrate and chemosynthetic fauna), 
and the end products (authigenic carbonates) of their degrada- 
tion are also known (Sassen et al, 2004). Anaerobic microbial 
oxidation of C2-C5 hydrocarbons has been inferred at the site 
from the enrichment in the S 13 C of the residual alkane pools 
(Sassen et al, 2004). In particular, geochemical measurements 
reveal a preferential degradation of propane, butane, and pentane 
(Sassen et al, 2004). 

To better understand the role of C2-C5 hydrocarbon degrading 
organisms in global geochemical cycles, we examined micro- 
bially mediated alkane consumption and SRRs, and the effect 
of alkane consumption on the inorganic carbon pool using sed- 
iments collected from a marine hydrocarbon-rich seep in the 
GoM. Specifically, we conducted a series of experiments in ex situ 
batch reactors to examine: ( 1 ) the rate at which microbial com- 
munities degrade C1-C4 alkanes; (2) the relationship between 
alkane degradation and sulfate reduction; (3) the degree to which 
microbially mediated alkane degradation influences the isotopic 
signatures of the alkanes and dissolved inorganic carbon (DIC) 
pools; and (4) how community composition of the GoM sedi- 
ments are affected by the addition of Q-C4 alkanes. This study 
advances our understanding by quantifying the potential rates of 
C1-C4 alkane consumption, sulfate reduction and the possible 
effects on the local carbon pool at a well-studied marine habitat. 
Furthermore, we describe the microbial phylotypes that are most 
abundant during active C1-C4 degradation. 

MATERIALS AND METHODS 
STUDY SITE AND SAMPLE COLLECTION 

Sediments were collected from the Garden Banks mud volcano 
site (GB425) in the northern GoM (27-33. 140 N, 92-32.437 W) 
at 597 m depth, during an expedition with the R/V Atlantis and 
DSVAlvin (Dive 4645) in November 2010. Intact sediment cores 
were recovered with polyvinylchloride core sleeves (20-30 cm 
height, 6.35 cm ID, 0.32 cm sleeve thickness). Sediment sampling 
sites were selected based on the presence of chemoautotrophic 
Beggiatoa mats overlying the sediments and the previous detec- 
tion of alkanes. Retrieved cores were sealed under Argon gas to 
limit gas exchange with the atmosphere and to prevent reoxi- 
dation of sulfide to sulfate, and refrigerated for transport to the 
laboratory. It is important to note that sulfide reoxidation is very 
rapid in these sediments, particularly in sediments hosting micro- 
bial mats (Bowles et al., 201 1). Therefore, it is not surprising that 
there is no observable sulfate gradient. Hydrogen concentrations 
at the study site were determined using a "reduction gas analyzer" 
as described by Orcutt et al. (2005). The in situ SRRs were mea- 
sured shipboard as described previously (Bowles et al., 201 1). The 
bottom water temperature at this study site was 8°C. 

BATCH REACTOR SET-UP AND SAMPLING 

Collected sediments were transferred to an anaerobic chamber 
with a 5% H 2 /75% N 2 /20% C0 2 atmosphere (Coy Laboratory 
Products, Grass Lake, MI). The sediment core used for this study 
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was stored for 3 months at 7°C. Sediments were thoroughly 
mixed, and were then diluted with an equal volume of anaerobic 
medium with 28 mM sodium sulfate and 2mM sodium sulfide 
(Widdel and Bak, 1992). No nitrate or nitrite was added to the 
medium. The resulting slurry was aliquoted into sterile 200 mL 
serum bottles and sealed with a butyl rubber stopper under strict 
anoxic conditions. The serum bottle headspace (100 mL) was 
flushed and then filled with a single gas at ~69 kPa (concentra- 
tions of 40-80 mM, Figure 1) of chemically pure (>99% purity) 
methane, ethane, propane, or butane (Airgas East, Waterford, 
CT, USA). The solubility at standard conditions (in water) of 
each gas is as follows: methane (0.9 mM), ethane (1.3 mM), 
propane (1 mM), and butane (0.8 mM) (webbook.nist.gov/) and 
exceeds the observed in situ concentrations of ~400, 20, 1, and 
1 u,M, respectively reported below. Excess gas pressure was used 
to overcome potential issues in isotopic data interpretation as 
reported for anaerobic propane oxidation previously (Quistad 
and Valentine, 201 1), and to avoid substrate limitation. The con- 
trol bottles were flushed and filled with chemically pure (>99% 
purity) nitrogen (Na) gas. An initial sample of each gas (except 
Na) was taken for isotopic analyses. The slurry was also sub- 
sampled and frozen at — 80° C for DNA extraction, sulfide and 
sulfate quantification. 

To test whether alkanes could be oxidized without concomi- 
tant sulfate reduction, we use the competitive inhibitor molybdate 
to inhibit sulfate reduction and monitored the subsequent alkane 
consumption rates. Briefly, 5 mL of sediment slurry was trans- 
ferred to 25 mL Balch tubes, in duplicate, per gas amendment, 
and sealed using a rubber butyl stopper inside an anaerobic cham- 
ber. These incubation volumes were used to maximize analyses 
given the limited sediment volumes. Due to the high sulfate 
concentration present in these incubations, the affect of changing 



the incubation volume should not be detrimental to sulfate 
reduction. Sodium molybdate was added to each tube to a final 
concentration of 28 mM (Orcutt et al., 2008). These tubes were 
then flushed and filled with the appropriate Q-C4 or N2 gas at 
~69 kPa. 

All the reactors were incubated at 7°C and the headspace was 
sampled every 15-days to monitor C1-C4 consumption. After 80 
days of incubation, the final gas concentrations were measured, 
and gas samples were archived in gastight Exetainers (Labco 
International, Houston, TX, USA) for natural abundance isotopic 
measurements. Samples were also withdrawn for DIC, DNA, sul- 
fide and sulfate measurement, and preserved by freezing at — 80° C 
in appropriate vials. All sub-samples were collected and measured 
in triplicate. 

GE0CHEMICAL MEASUREMENTS 

C1-C4 alkanes were quantified from the headspace by subsam- 
pling a 50 (AL aliquot and analyzing alkane concentrations on a 
gas chromatograph (Hewlett Packard 5890 Series II) equipped 
with a flame ionization detector and a packed column (RestekRt- 
XL). Chemically pure alkanes (>99% purity) (Airgas East, 
Waterford, CT, USA) were used to generate standard curves. To 
account for potential alkane leakage from the bottles, we set up 
sediment-free controls and monitored changes in alkane con- 
centration over time. Sediment-free controls showed ~4 ± 1% 
variation in gas measurements, which represents both the ana- 
lytical resolution of our measurements and/or modest loss of 
gas due to leakage or sorption into the stoppers, and is well 
below the rates of loss observed in the biological treatments 
(see below). Sulfate concentrations were determined using the 
QuantiChrom™ Sulfate Assay Kit (BioAssay Systems, Hayward, 
CA, USA). Sulfide concentrations were measured using a colori- 
metric assay based on the Cline method (Cline, 1969). Nitrate was 
measured using the resorcinol method as described previously 
(Zhang and Fischer, 2006). Nitrite was measured as previously 
described (Pai et al, 1990). 

MICROBIAL SULFATE REDUCTION RATE MEASUREMENTS 

SRRs were measured using a previously described radiotracer 
method (Fossing and Jorgensen, 1989). Briefly, the slurry incu- 
bations were opened under anaerobic conditions and 5 mL sub- 
samples of the enrichments were transferred to a Balch tubes. 
The tubes were sealed and pressurized as described previously. 
The Balch tubes were then amended with ca. 10 u,L of Na 35 SO^~ 
(2 (jiCi) and incubated for 1 day. Following incubation the slurries 
were shaken and 1 mL of slurry was drawn by syringe into 5 mL 
of 20% zinc acetate and shaken, effectively trapping H^Sfaq) as 
ZnS( s ). The ZnS solution was placed into a 15 mL Falcon tube, 
and washed three times with a 3% NaCl solution to remove any 
residual 35 SO^. Sulfide was extracted using the hot chromium 
reduction method (Fossing and Jorgensen, 1989), ultimately trap- 
ping sulfide in 20% Zn-acetate. The activity of 35 S was determined 
by liquid scintillation and SRR were calculated after Fossing and 
Jorgensen (1989) using Equation 1. 

T ? 1 1-06 a , , 

srr= sor * — * 1 

L 4 J t (A + a) 



Methane 



B Ethane 




C Propane 



D Butane 





FIGURE 1 I Alkane consumption as a function of time by the Gulf of 
Mexico (GoM) site GB425 (dive 4645 and core 22) incubations. The 

alkane concentration in the headspace of the sediment incubations was 
measured using a Hewlett Packard 5890 Series II gas chromatograph 
equipped with a flame ionization detector, on a Restek Rt-XL Sulfur packed 
column. (A) methane; (B) ethane; (C) propane; and (D) butane. Rates of 
alkane consumption were calculated using all available time points based 
on a linear regression. 
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Where [SO4 - ] is the concentration (nmol mL -1 ) of sulfate incu- 
bation, a is the activity (dpm) of the trapped sulfide, 1.06 is the 
fractionation factor between the sulfide and sulfate pools, A is the 
activity of the sulfate pool, and f is the incubation time (days). 
The rates are presented in units of nmol S mL~ 1 day - 1 . 

ISOTOPE ANALYSIS 

All isotopic analyses were performed at the Stable Isotope Facility 
at University of California, Davis using a method modified from a 
previous publication (Atekwana and Krishnamurthy, 1998). For 
DIC measurements, 1 mL filtered (0.2 |im) water samples were 
collected and injected into evacuated 12 mL septum capped vials 
(Exetainers, Labco, Houston, TX, USA) containing 1 mL 85% 
phosphoric acid. The evolved CO2 was purged from vials through 
a double-needle sampler into a helium carrier stream (20 mL 
min -1 ). For high concentration samples, gases were sampled by 
a six-port rotary valve (Valco, Houston, TX, USA) with a 100 |xL 
loop programmed to switch at the maximum CO2 concentration 
in the helium carrier. For low concentration samples, the entire 
CO2 content was frozen in a trapping loop then released to the GC 
column. The CO2 was passed to the IRMS through a Poroplot Q 
GC column (25 m x 0.32 mm ID, 45°C, 2.5 mL/min). A reference 
CO2 peak was used to calculate provisional delta values of the 
sample CO2 peak. Final 8 13 C values were obtained after adjusting 
the provisional values such that correct 8 13 C values for laboratory 
standards were obtained. Two laboratory standards were analyzed 
every 10 samples. The laboratory standards are lithium carbon- 
ate dissolved in degassed, deionized water, and a deep seawater 
reference material (both calibrated against NIST 8545). 

For isotopic analyses of the C1-C4 gases, a ThermoScientific 
PreCon concentration system interfaced to a ThermoScientific 
Delta V Plus isotope ratio mass spectrometer (Thermo- 
Scientific, Bremen, DE) was used as described previously (Yarnes, 
2013). Gas samples were purged from Exetainers through a 
double-needle sampler into a helium carrier stream (20 mL/min), 
which is passed through a H2O/CO2 scrubber [Mg(C104)2, 
Ascarite] and a cold trap cooled by liquid N2. The gas was sep- 
arated from residual gases by a Rt-Q-BOND GC column (30 m x 
0.32 mm x 10 u,m, 30°C, 1.5 mL/min). After the gas eluted from 
the separation column, it was either oxidized to CO2 by reaction 
with nickel oxide at 1000°C (8 13 C), or pyrolyzed in an empty alu- 
mina tube heated to 1350°C (8 2 H) and subsequently transferred 
to the IRMS. A pure reference gas (CO2 or H2) was used to cal- 
culate provisional delta values of the sample peak. Final 8-values 
are obtained after adjusting the provisional values for changes 
in linearity and instrumental drift such that correct 8-values for 
laboratory standards were obtained. Laboratory standards were 
commercially prepared gases diluted in helium or air and were 
calibrated against NIST 8559, 8560, and 8561. 

DNA EXTRACTION, MASSIVELY PARALLEL SEQUENCING. AND 
PHYL0GENETIC ANALYSES 

Sediment was subsampled under anoxic conditions for To and Tf 
for nucleic acid extractions. These samples were flash frozen in 
liquid N2 and stored at — 80° C until use. DNA was extracted using 
the PowerSoil® DNA Isolation Kit (MO BIO Laboratories, Inc., 
Carlsbad, CA, USA) as per the manufacturer's guidelines. The 



extracted DNA was subjected to massively parallel sequencing of 
the 16S ribosomal RNA (rRNA) gene using Roche 454 Titanium™ 
chemistry and the primer pairs 27F/519R and 340F/806R for the 
bacterial V1-V3 and archaeal V3-V4 regions, respectively (Dowd 
et al., 2008; Acosta-Martinez et al., 2010). The resulting sequences 
were analyzed as previously described, and denoised using the 
QIIME pipeline (Adams et al, 2013). Phylogenetic analysis was 
performed as previously described using FastTree(2.1.7) for tree 
generation with 25 representative sequences (Adams et al, 2013). 
All sequences generated in this study are deposited with NCBI 
(accession #SRP032824). 

QUANTITATIVE-PCR 

Quantitative PCR (qPCR) was used to determine the abun- 
dance of bacterial and archaeal 16S rRNA, dsrA, aprA, and mcrA 
genes. In addition, qPCR was used to enumerate the abundance 
of sulfate-reducing prokaryotes by amplifying the adenosine 5'- 
phosphosulfate [APS] reductase (aprA) gene with primers spe- 
cific to sulfate-reducing bacteria and archaea (Christophersen 
et al, 2011). Primers specific to the bacterial dissimilatory sul- 
fite reductase (dsrA) gene were used to quantify members of 
sulfate-reducing bacteria (Kondo et al, 2004). We refer to all 
sulfate-reducing microbes as sulfate-reducing prokaryotes or 
SRP throughout. Methanogenic archaea were quantified using 
mcrA primers directed specifically toward the methanogenic 
methyl Coenzyme M reductase encoding gene (Luton et al, 2002; 
Ver Eecke et al., 2012). Standard curves were constructed by 
serial dilution of linearized plasmids containing the target gene 
(Table 1). Quantification was performed in triplicate with the 
Stratagene MX3005p qPCR System (Agilent Technologies) using 
the Perfecta SYBR FastMix with low ROX (20 u.L reactions, 
Quanta Biosciences, Gaitherburg, MD) and specific primers and 
annealing temperatures (Table 1). The temperature program for 
all assays was 94°C for lOmin, 35 cycles of 94° C for 1 min, the 
annealing temperature for 1 min (Table 1), extension at 72°C for 
30 s, and fluorescence read after 10 s at 80° C. Following amplifi- 
cation, dissociation curves were determined across a temperature 
range of 55-95°C. C f -values for each well were calculated using 
the manufacturer's software. 

RESULTS 

GE0CHEMICAL CHARACTERISTICS OF SITE GB425 

Though the gross geochemistry of this site has been previously 
described (Joye et al., 2009), here we present the alkane concentra- 
tions and other geochemical attributes of the specific sediments 
used in these studies (Table 2). DIC ranged between 4 and 6 mM, 
while the observed dissolved organic carbon (DOC) is about 
1-3 mM through the sediment depths surveyed. Nitrate and 
nitrite concentration was 5-40 u,M in the upper layers of the sed- 
iment. Sulfate, the dominant oxidant, was replete throughout the 
sediment profile (24-36 mM) and was higher than typical seawa- 
ter values (28 mM) (Canfield and Farquhar, 2009). n-alkanes were 
observed only between 9 and 15 cm sediment depth. Between the 
depth ranges of 9-12 and 12-15 cm ethane was observed at 17.22- 
22.33 ^M, propane at 1.45-0.75 (xM, butane at 0.74-0.35 [iM. 
Pentane was not observed. Methane concentrations peaked at 
425.05 [lM, at ~15-18 cm sediment depth (Table 3). 
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Table 1 | Primers and Conditions for quantitative PCR assays. 



Process 


Target gene 


Forward primer (nM) 


Reverse primer (nM) 


Positive control 


Annealing 
temp in °C 
(References) 


Sulfate reduction 


Dissimilatory 

sulfite 

reductase 


DSR1-F+W0) 

ATCGGNCARGCNTTYCCNTT 


DSR-R (600) 

GTGGMRCCGTGCAKRTTGG 


Desulfovibrio 
vulgaris 


58° C (Kondo 
et al., 2004) 


Sulfate reduction 


Adenosine 5'- 

phosphosulfate 

reductase 


aps3F (400) 

TG G C AG ATCATG WTYAAYG G 


aps2R (400) 

GCGCCGTAACCRTCYTTRAA 


Desulfovibrio 
vulgaris 


55° C 

(Christophersen 
etal., 2011) 


Methanogenesis 


Methyl CoM 
reductase 


qmcrAF-alt (150) GAR GAC 
CACTTYGGH GGTTC 


MLR (200) 

TTC ATTG C RTAGTTWG G RTAGTT 


Methanosarcina 
acetivorans, 
Methanococcus 
jannaschii 


59°C (Luton 
et al., 2002; 
Ver Eecke 
et al., zo\Z) 


Bacteria 


16S rRNA 


Bact1369F (1000) GTT GGG 
GCC RCC WCK KCK NAC 


Prok1541R (1000) 
CGGTGAATATGCCCCTGC 


Arcobacter 
nitrofigulis 


59° C (Suzuki 
etal., 2001) 


Archaea 


16S rRNA 


Arch349F (500) 
GYGCASCAGKCGMGAAW 


Arch806R (500) 

G G ACTAC VS G G GTATCTAAT 


Ferroplasma 
acidarmonas 
Fer1 


54° C (Takai 
and Horikoshi, 
2000) 



Table 2 | Geochemical data from site GB425, from which sediments were collected for these analyses in November 2010 (27°33.1887N, 
93°32.4449W). 

Depth DIC d 13 C- Hydrogen DOC Sulfate Sulfide Nitrate Methane pH Sulfate Anaerobic 

range (mM) DIC (nM) (mM) (mM) (mM) and Nitrite (nM) reduction rate methane oxidation 

(cm) (%c) (mM) (nmol mL" 1 day -1 ) (nmol mL" 1 day -1 ) 



Overlying water 


n.s. 


n.s. 


n.s. 


958 


n.s. 


n.s. 


23.1 


n.s. 


7.5 


n.a. 


n.a. 


0-3 


3.9 


-10.0 


37.4 


2477 


24.2 


0.6 


36.2 


122.3 


7.5 


86 ±5 


0.8 ± 0.3 


3-6 


4.4 


-10.5 


242.4 


1203 


32.3 


0.1 


5.2 


63.0 


7.6 


344 ± 216 


1.6 ± 0.3 


6-9 


4.2 


-12.6 


21 


1702 


32.3 


0.6 


b.d.l. 


378.1 


7.9 


182 


6.8 ± 0.1 


9-12 


5.1 


-19.5 


29.8 


1.S. 


27.9 


0.2 


b.d.l. 


371.9 


i.v. 


135 ±22 


5.0 ± 2.5 


12-15 


4.7 


-15.2 


51.9 


1.S. 


n.s. 


0.8 


b.d.l. 


413.1 


i.v. 


n.s. 


n.s. 


15-18 


4.9 


-20.1 


31.4 


1.S. 


29.9 


0.6 


b.d.l. 


425.1 


8.0 


596 ± 8 


30.3 ± 7.6 


18-21* 


5.8 


-28.6 


38.9 


2618 


29.1 


n.s. 


40.6 


200.2 


i.v. 


185 ±216 


7.8 ± 0.5 


21-24 


6.0 


-23.7 


39.6 


2982 


25.0 


0.1 


b.d.l. 


357.7 


i.v. 


743 ± 340 


18.2 ± 10.0 


24-27 


5.3 


-23.3 


50.7 


1420 


35.8 


1 


b.d.l. 


288.8 


i.v. 


287 ± 120 


18.9 ±0.4 


27-30 


5.3 


-24.4 


31.5 


1.S. 


n.s. 


1.2 


b.d.l. 


203.2 


i.v. 


n.s. 


n.s. 



DIC, Dissolved inorganic carbon; DOC, Dissolved organic carbon, n.s., No sample; l.s., Lost sample; i.v., Insufficient volume; b.d.l., Below detection limit; n.a., Not 
applicable. Rate measurements are mean ± standard error (n = 2). 'Data at this depth range appears unreliable. 



Ci-C 4 ALKANE OXIDATION OCCURS IN BATCH REACTORS 

n-Alkane consumption began within the first 15 days of the 
80-day incubations for C1-C3 gases (defined as >10% con- 
sumption compared to To) (Figure 1). C4 consumption was 
only measurable after ~45 days of incubation. The high- 
est % consumption was observed for propane (73 ± 13%) 
followed by butane (45 ± 5%), ethane (31 ± 6%), and 
methane (25 ± 6%). The highest rate of consumption was 
observed for propane (354 ± 37 nmol mL~' day -1 ) fol- 
lowed by methane (263 ± 68 nmol cm -3 day -1 ), ethane (168 
± 5 nmol cm -3 day -1 ), and butane (125 ± 16 nmol cm -3 
day -1 ) (Figure 2, Table 4). Along with alkane oxidation we also 



observed a decline in sulfate concentrations and a concomitant 
increase in sulfide concentrations (Table 5). Importantly, the 
sulfide concentrations were below those observed to be 
inhibitory (16.1 mM) for sulfate-reducing bacteria (Reis et al., 
1992). 

SULFATE REDUCTION IS COUPLED TO C 1 -C 4 ALKANE OXIDATION 

The addition of each Q-C4 gas increased the SRR over the N2 
control treatment by at least 2-fold (Figure 2). The rates reported 
from these incubations are comparable to previous reports from 
GoM non-seep porewaters and sediments (Arvidson et al., 2004; 
Joye et al., 2004) but lower than those measured shipboard on 
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freshly collected samples (Table 2). While it is impractical to iden- 
tify the precise cause of this discrepancy, there are a few likely 
factors that could have contributed to these differences, includ- 
ing (A) natural heterogeneity in the geochemistry and microbial 
community composition and activity; (B) the process of sediment 
homogenization prior to incubations, which does not represent 



Table 3 | C1-C5 alkane concentrations in sediments at site GB425, 
from which sediments were collected for this study in November 
2010 (27°33.1887N, 93°32.4449W). 



Depth range 


Methane 


Ethane 


Propane 


n-butane 


Pentane 


(cm) 


(lilvl) 


<M.M) 


<M.M) 


(|iM) 


(^M) 


9-12 


371.94 


17.22 


1.45 


0.74 


0.00 


12-15 


413.12 


22.33 


0.75 


0.35 


0.00 




ra 



400- 



E 300- 




<? <? <? c& 

^ 



FIGURE 2 | Potential sulfate reduction rates (SRR) were measured 
using the 35 SO^~ radiotracer method (Fossing and Jorgensen, 1989) 
and consumption rates for C1-C4 alkanes by alkane amended slurries 
of GoM site GB425 sediments. The SRR assays were performed for 24 h. 
Values represent average ± standard deviation of triplicate measurements 
of duplicate incubations. Alkane consumption rates were calculated from a 
linear regression as in Figure 1. Rates of sulfate reduction were calculated 
as described in the methods. 



maximal or minimal rates; and (C) changes in microbial compo- 
sition and activity during the 3 months of storage. Comparison 
of SRR in the one-day incubations to C1-C4 alkane oxidation 
rates (Table 4) shows that the addition of methane or any of the 
four alkanes stimulates SRR over the N2 control treatment. From 
the predicted reaction stoichiometry (Table 6) both the methane 
and ethane oxidation rates correspond closely with the observed 
SRR. In contrast substantially higher levels of propane and butane 
oxidation were observed than can be supported by sulfate reduc- 
tion alone. Incubation with molybdate inhibited the oxidation of 
C1-C4 gases by ~90-97% (Table 4), consistent with the direct 
involvement of sulfate-reducing prokaryotes in alkane oxidation. 



CARBON IS0T0PIC SIGNATURE AND ANALYSES 

The 8 13 C signature of the methane in the headspace did not 
appreciably change over the course of the incubation period 
(Table SI). This contrasts with the isotopic signatures of the other 
alkanes. As mentioned, ethane concentration decreased to about 
10 mM over the course of the incubation, but the isotopic change 
in the pool was not significant. The incubations with propane 
showed the largest decrease in concentration (~30mM) over the 
80-day incubation period. Over this time, the propane pool was 
enriched by 4.4%o. Finally, incubations with butane resulted in a 
decrease in the pool size of ~ 10 mM ( 1 mmol) and an enrichment 
in the residual butane pool of 4.5%o. 



Table 5 | Sulfate and sulfide concentrations measured in the initial 
sediment slurry and at the final time-point. 

Sulfate concentration Sulfide concentration 
(mM) (mM) 



Sediment slurry (Initial) 31 .6 ±1.2 

N 2 control (Final) 26.4 ±1.8 

Methane (Final) 26.1 ±6.1 

Ethane (Final) 22.1 ±4.0 

Propane (Final) 15.6 ±4.1 

Butane (Final) 13.3 ±1.2 



2.3 ±0.2 
8.2 ± 0.1 
9.1 ± 0.4 

10.4 ± 1.1 
9.5 ± 1.4 

15.4 ±3.2 



Table 4 | Comparing rate of alkane oxidation and sulfate reduction, and the effect of molybdate on alkane oxidation. 





Rate of alkane 
consumption 
(nmol mL" 1 
day" 1 ) 


Rate of alkane 

consumption 

with 

molybdate 
addition (nmol 
mL" 1 day 1 ) (% 
inhibition) 


Total sulfate 
reduction rate 
(nmol mL -1 
day" 1 ) 


Observed 
ratio using 
total SRR 
alkane:SO|" 


Corrected sulfate 
reduction rates 
above the nitrogen 
control (nmol mL -1 
day" 1 ) 


Observed ratio 
using 

corrected SRR 
alkane:SO|" 


Predicted 
ratio of 
alkane: SO|" 
from Table 6 


Methane 


263 ± 68 


29 ± 1 (89) 


297 ± 33 


0.9 


203 ± 53 


1.3 


1 


Ethane 


168 ± 5 


12 ±4 (93) 


330 ± 121 


0.51 


236 ± 127 


0.7 


0.57 


Propane 


354 ± 37 


10 ±5 (97) 


246 ± 44 


1.4 


152 ± 60 


2.3 


0.4 


Butane 


125 ± 16 


14 ±2 (89) 


220 ± 38 


0.57 


125 ± 56 


1 


0.34 



Rates of alkane consumption were calculated using all available time points based on a linear regression. Rates of sulfate reduction calculated as described in the 
methods. 
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Table 6 | Gibbs free energy of the anaerobic oxidation of acetate, methane, and alkanes using sulfate as an electron acceptor (conditions 
shown are at standard temperature and pressure). 





Sulfate reduction process 


Reaction 


zivj |KJ/niol Ow^ J 


Csrbon sourc&i SO^ 




1 


Heterotrophic (acetate) 


SOj" + CH3COO- -► 2HCO3 + HS- 


-47.7 


1:1 


2:1 


2 


Heterotrophic (methane) 


SOj" + CH 4 ->• HCO3 + HS" + H 2 0 


-33 


1:1 


1:1 


3 


Heterotrophic (ethane) 


14SO;;" + 8C 2 H 6 ->• 14HS" + I6HCO3- + 8H 2 0 + 2H+ 


-39.81 


8:14 


16:14 


4 


Heterotrophic (propane) 


5SO^- + 2C 3 H 8 -»• 6HCO3- + 5HS" + H+ + 2H 2 0 


-33.06 


2:5 


6:5 


5 


Heterotrophic (butane) 


26SO4" + 9C 4 H 10 + 4H 2 0 -»• 36HCO- + 36H+ + 26HS" 


-14 


9:26 


18:13 


6 


Autotrophic (with H2) 


SOj" + 2HCO3 + 8H 2 + 2H+ -+CH 3 COCr + HS" + 8H 2 0 


-336.5 


2:1 


2:1 



Autotrophic sulfate reduction, in which hydrogen is used to reduce inorganic carbon, is shown for reference. 

*AG° Values reported are those calculated under standard conditions of 7 M concentrations for soluble reactants, 7 atmosphere pressure for gases, 298. 15 K 
temperature at pH 7.0 and are calculated using values from the CHC Handbook for Chemistry and Physics (http://www.hbcpnetbase.com/). 
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FIGURE 3 | Relative abundance (percentage) of Proteobacteria and 
Euryarchaeota determined from massively parallel pyrosequencing of 
DIMA extracted from batch reactor sediments incubated with methane, 
ethane, propane, butane, and nitrogen and pre-incubation [To) 
sediments. Top (A) and bottom (B) panels show the taxonomic breakdown 
of sequences at the class and order level, respectively. Sequences sharing 
97% nucleotide sequence identity are defined as operational taxonomic 
units (OTUs). 



MICROBIAL COMMUNITY ANALYSES 
454 pyrotag sequencing 

A total of 11,725, 17,003, 12,529, 16,208, and 18,015 bacterial 
sequences were analyzed from sediments incubated with N2, 
methane, ethane, propane, and butane, respectively, and 12,944 
bacterial sequences from the Tq sediment. There were shifts 
between the Proteobacterial communities of the alkane batch 
reactors in comparison to the control and To sediment com- 
munity (Figure 3A). Among sequences allied to known sulfate - 
reducing Deltaproteobacteria, there was an increase from the 
To sequences (~20%) in the N2, methane, ethane, propane, 
and butane sequence libraries (~23, 32, 23, 33, and 55%, 
respectively) (Figure 3A). In turn, there was a decrease in the 
representation of Gammaproteobacteria in the N2, methane, 
ethane, propane, and butane sequence libraries (~37, 17, 30, 
26, and 12%, respectively) from the Tq sequences (~53%). 16S 
rRNA gene phylogeny revealed that the ethane reactors har- 
bored a putative SRP community that was distinct from the 
propane and butane reactors (Figure 4). These sequences com- 
prised the majority (90-95%) of the Deltaproteobacterial com- 
munity (Figure 4). In the ethane reactor community, the most 
closely related Deltaproteobacterial 16S rRNA gene sequences 
(95-99% nucleotide sequence identity) included strain BuS5 
(accession no. EF077225), the enrichment culture "Butanel2- 
GMe" (accession no. EF077226), and other SRP clones from 
sediments retrieved from the GoM (clone GoM_DSSGM3_28, 
accession no. FR872064; clone GoM_DSSGM3_19, accession no. 
FR872059; and clone GoM161_Bac9, accession no. AM745163) 
(Kniemeyer et al., 2007; Orcutt et al, 2010; Kleindienst et al., 
2012). In contrast, SRP sequences in the propane and butane 
batch reactor communities were most closely allied to uncul- 
tured Deltaproteobacteria clones from propane- and butane- 
oxidizing enrichments of hydrocarbon seep sediments from the 
GoM (Propanel2-GMe clone 230, accession no. FR823371) and 
Hydrate Ridge (Butanel2-HR clone 302, accession no. FR823375 
and Butanel2-HR clone 342, accession no. FR823377) (Jaekel 
etal., 2013). 

A total of 18,667, 10,291, 18,545, 12,462, 9743, and 13,233 
archaeal sequences were also analyzed from the N2, methane, 
ethane, propane, and butane batch reactors and To sediments, 
respectively. There were notable shifts in the sequences allied to 
the class Methanomicrobia from the initial sediment community 



and across the different alkane batch incubations (Figure 3B). 
Over 58% of sequences were allied to Methanomicrobia in To 
sediments, increasing to comprise ~87 and 94% of methane 
and butane sequences. Within the Methanomicrobia, there were 
also notable changes in sequences identified as phylotypes that 
mediate AOM. For the putative methane-oxidizing communi- 
ties, ANME-1 comprised ~40% of the Methanomicrobia in the 
incubation with methane, but less than 5% of sequences were 
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FIGURE 4 | Phylogenetic affiliation of 16S rRNA Deltaproteobacterial 
gene sequences retrieved from Gulf of Mexico batch reactor 
sediments. A total of 25 representative sequences from Gulf of Mexico 
sediments incubated in batch reactors with ethane (GOM_ETHANE), 
propane (GOM_PROPANE), and butane (GOM_BUTANE) are shown in bold. 
The phylogenetic tree was generated by maximum likelihood with FastTree 
Version 2.1.3. Local support values shown are based on the 
Shimodaira-Hasegawa (SH) test with 1000 resamples. Only values >80% 
are shown on the branches as black circles. The 16S rRNA sequence of 
Archaeoglobus profundus DSM 5631 (NR_074522) was used as an 
outgroup. Scale = 0.1 substitutions per site. 



allied to ANME-1 in the To, N2, ethane, propane, and butane 
sediments (Figure SI). 

Quantitative PCR 

qPCR using specific primers for 16S rRNA showed that bacte- 
rial 16S rRNA gene abundance was two orders of magnitude 
higher than archaeal 16S rRNA gene abundance at the start of 
the incubation (Figures 5IA,B)- Bacterial abundance was only 
slightly elevated (less than an order of magnitude) over the Tq 
assessment in all treatments at the end of the incubation period 
with the greatest increase in population observed in the N2 and 



CH4 amendments (~3-fold increase). Addition of alkanes also 
stimulated bacterial population growth of about 2-fold over the 
initial population estimates. N2 and ethane amendments resulted 
in a 3-fold increase in archaeal populations while propane and 
butane yielded a 1.5-fold increase. These differences are con- 
sistent among treatments. However, 40% of bacterial genomes 
contain 1-2 copies of rRNA genes, though microbial genomes 
with as many as 15 copies have been reported (Acinas et al., 
2004). Moreover, archaeal genomes are known to harbor between 
1 and 5 rRNA gene copies per genome (Acinas et al, 2004). Thus, 
given these differences, as well as environmental heterogeneity 
and other factors, the differences presented here likely reflect rel- 
ative changes in proportion, but the significance of these changes 
among treatments remains unconstrained. 

Estimates of aprA gene abundance, a marker for SRP, reveal the 
highest abundances at the initiation of the incubation and under 
the N2 amendments. aprA gene abundance after incubation with 
methane or the tested alkanes show a decrease in abundance, per- 
haps indicating a shift in community due to incubation effects 
that is consistent with the decrease in bacterial 16S rRNA gene 
abundance over the same treatments. Notably, of the alkane addi- 
tions, propane maintained the largest SRP population followed by 
ethane, methane and butane treatments (Figure 5IIA). 

The 16S rRNA gene phylogeny of SRP is diverse and difficult to 
capture with specific primers. Accordingly, we employed primers 
targeting the gene encoding for dissimilatory sulfite reductase 
(subunit A) with primers that target both Gram-positive and 
Gram-negative bacterial species of SRP (Kondo et al., 2004). 
With the exception of the N2 treatment, the dsrA gene abun- 
dance is similar across all treatments. The dsrA gene abundance 
in the N2 treatment is about 3-fold higher than observed in other 
treatments (Figure 5IIB). 

Archaeal mcrA gene abundance was highest, and had the great- 
est standard deviation at the initial sampling. The lowest observed 
mcrA gene abundance occurred in the methane amended samples 
and concurs with a similar observed decrease in the total archaeal 
16S rRNA gene abundance described above. This suggests that the 
addition of methane adversely affected the natural methanogen 
population over the course of the incubation (Figure 5IIIA). 

DISCUSSION 

The anaerobic microbial degradation of short-chain alkanes has 
recently gained attention because microbes mediating these pro- 
cesses may compete for the oxidant pool (sulfate), potentially 
influencing the rates of AOM (Kniemeyer et al, 2007; Savage 
et al, 2010; Adams et al., 2013; Jaekel et al, 2013). The data 
herein reveal that C1-C4 alkane consumption — including anaer- 
obic ethane oxidation — stimulated sulfate reduction. To assess the 
significance of the determined potential rates of sulfate reduction 
and w-alkane consumption, two points must be addressed: (1) the 
intrinsic sulfate reduction activity in the GoM sediments used, 
and (2) the concentrations of the substrates used relative to those 
measured in situ. First, bulk geochemical analyses show that the 
GoM sediments are rich in organic matter and hydrogen that can 
support the growth of heterotrophic and autotrophic SRP. At this 
study site, the DOC concentration was ~l-3 mM and hydrogen 
was in the nM range. It was therefore critical to account for sulfate 
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FIGURE 5 | Abundance of microbes determined using quantitative 
PCR. Panel (I) represents the 16S rRNA abundances for A: Bacteria 
and B: Archaea. Panel (II) represents the abundance of 



sulfate-reducing bacteria as determined using A: aprA and B: dsrA. 
Panel (III) represents the abundance of methanogens as determined 
using A: mcrA. 



reduction attributable to endogenous electron donors, which we 
did by maintaining the native sediments under a N2 atmosphere. 
Not surprisingly, these incubations exhibit intrinsic sulfate reduc- 
tion activity at nearly lOOnmol SO^ mL _1 day -1 . Importantly, 
the addition of n-alkanes increased this baseline sulfate reduction. 
Second, sulfate concentrations used in the ex situ incubations cor- 
respond with those measured at various depth ranges at this GoM 
marine seep (Table 2) and were not limiting (in bulk geochem- 
istry; Table 5) over the course of the incubation. Given that these 
experiments were conducted at conditions that might favor anaer- 
obic alkane oxidation, e.g., an abundance of one alkane in the gas 
phase, and a media replete with sulfate, these data provide insight 
into the relationship between alkane oxidation and sulfate reduc- 
tion, and represent "potential" rates of n-alkane consumption 
linked to sulfate reduction. 

With the aforementioned points in mind, we address the link- 
age between alkane consumption and sulfate reduction from a 
few perspectives. First, based on the stoichiometry of each alkane 
oxidation-sulfate reduction pathway (Table 6) we estimate the 
contribution of each oxidation pathway to SRR. Second, we use 
the change in the isotopic signatures of the alkanes and the 
DIC pool to estimate carbon exchange between the alkane and 
DIC pools, and compare this carbon mobility with the alkane- 
oxidation rates (Table S2). Furthermore, we examine the commu- 
nity composition within the enrichments in an effort to elucidate 



community members potentially responsible for alkane oxidation 
and sulfate reduction (Figures 3-5). 

From the stoichiometry of the reaction pathways (Table 6), 
which assume the alkanes are completely oxidized to HCO3 , and 
that no alkane-derived carbon is assimilated, the linkage between 
alkane oxidation and sulfate reduction can be estimated. The 
methane incubation showed a potential net consumption rate of 
263nmolCmL _1 day -1 (Table 4). This coincides with a potential 
SRR of 297nmol S mL~' day -1 resulting in a carbon to sul- 
fate ratio of (0.9), consistent with the stoichiometric prediction 
of C:S of 1 . These estimates assume that all the sulfate reduction 
observed in the incubations is a result of methane oxidation, how- 
ever, the N2 control treatments indicate a potential intrinsic SRR 
of 94 nmol mL -1 day -1 . If we assume changes in the community 
activities between the methane addition and control treatments 
are negligible, aside from the oxidation of the methane, then 
we can use the SRR of the control treatment as a background 
SRR. Correcting for the intrinsic sulfate reduction results in an 
apparent SRR of 203 nmol mL _1 day -1 and an excess methane 
consumption (a C:S of 1.3, Table 4). 

One way in which more carbon may be consumed than pre- 
dicted from stoichiometry is if the alkane is not completely 
oxidized to DIC. The shifts in the isotopic composition of the DIC 
and alkane pools can be used to constrain how much carbon has 
moved from the alkane pool into the DIC pool and thus establish 
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if another carbon sink may be important. The moles of carbon 
from the alkane pool needed to shift the DIC pool from its initial 
to final composition can be described by: 

o 13 C-DIC T / * [DICto + Alk] * V = o 13 C-DIC T o 
* [DICto] * V + o 13 C-Alk T0 * [Alk] * V (2) 

where 8 13 C-DIC and S 13 C-Alk represents the isotopic signature 
of the pool, [DICro] is the concentration of the initial DIC pool 
(molCL -1 ), [Alk] is the alkane carbon oxidized (molCL -1 ), and 
V is the incubation volume (100 mL). Initial [DIC7-0] and S 13 C- 
DlCro were calculated as a 1:1 (v:v) mixture of measured values 
from the media and average values from the pore water (Table 2) 
which comprised the slurry incubations. Using these values the 
amount of carbon transferred from the initial alkane pool can 
be estimated. For the methane incubations, this calculation sug- 
gests that 2.2 mmols C are transferred from the methane to the 
DIC pool over the course of the incubation, indicating a rate of 
methane consumption of 275 nmol mL -1 day -1 . This is reason- 
ably consistent with the rate estimates derived from the change in 
methane concentration over time. 

When ethane oxidation is coupled to sulfate reduction, it 
results in a carbon to sulfate ratio of 8:7 (or 0.57, Table 6). 
The reported oxidation rate of 168 nmol ethane mL -1 day -1 
(or 336 nmol C mL -1 day -1 ) would lead to the consumption of 
294 nmol SO^ - mL - 1 day - 1 . This oxidation rate accounts for the 
majority (89%) of the total estimated sulfate reduction (330 nmol 
SO4 - mL -1 day -1 ) during the incubation. When SRR are cor- 
rected for the potential intrinsic rates, it results in a corrected rate 
of 236 nmol SO^ - mL -1 day -1 . The observed ratio of carbon 
oxidized to sulfate reduced is 0.7, about 25% more carbon than 
expected (Table 4). Unfortunately, the analytical resolution of the 
S 13 C-ethane precludes estimation of carbon movement among 
pools (Table SI). 

In the case of propane, the estimated SRR increased 2.6- 
fold over the N2 control, with an observed oxidation rate of 
354 nmol propane mU 1 day -1 (or 1062 nmol C mU 1 day -1 ). 
Given the stoichiometric relationship of 6:5 (Table 6), this con- 
sumption corresponds to 885 nmol SO4 - mU 1 day -1 , a rate 
3.6-fold higher than the measured potential SRR (246 nmol S 
mL -1 day -1 ) (Table 4). However, examining the change in the 
DIC pool, we calculate that 3.9 mmols C (1.3 mmols of propane) 
moved from the propane pool to the DIC pool (Table SI). This 
is equivalent to an oxidation rate of 487 nmol C mL -1 day -1 , 
accounting for 46% of the total loss of propane. Thus, there must 
be another sink for propane (discussed below), implicating the 
presence of another oxidant, another source of light carbon to the 
DIC, or both. 

Similar to the propane treatment, the butane addition also 
resulted in higher SRR (2.3-fold) compared to the control treat- 
ment. Over the course of the experiment, butane was consumed 
at a rate of 125 nmol butane mL -1 day -1 (or 500 nmol C mL -1 
day -1 ) (Table 4). The corresponding total (220 nmol S mL -1 
day -1 ) and corrected SRR (125 nmol S mL -1 day -1 ) can only 
account for about one-third of the butane consumption (Table 4). 
However, unlike propane, the majority (85%) of the butane 



carbon can account for the change in the DIC pool (Table SI). 
Again it is possible that another sink may exist but identifying 
such a sink is beyond the scope of these data. 

The data above underscore that the rate of oxidation of 
ethane, propane, and butane cannot be explained solely by the 
estimated rates of sulfate reduction. Similar observations have 
been reported by Quistad et al., who noted that the propane 
loss they observed might be accounted for by abiotic processes 
such as leakage and dissolution, partial degradation to alco- 
hols or acids, and/or inaccuracy in measurements (Quistad and 
Valentine, 2011). In our study, however, the observed higher rate 
of n-alkane oxidation may be best explained by (1) utilization of 
oxidants other than sulfate (e.g., NO^, which is present at site 
GB425 though not measurable in our reactors; Table 2), (2) errors 
in the estimation of either the oxidation or reduction rates or iso- 
topic assays due to systemic errors, (3) changes in the microbial 
community or activity of the community over the course of the 
incubation that were not observable with the sampling design, or 
(4) the precipitation of (authigenic) carbonate in the batch reac- 
tors as has been noted to occur in GoM sediments (Sassen et al., 
2004). We address each of these possibilities in detail below. 

Nitrate and/or nitrite represent potential alternative oxidants 
for alkane oxidation, and are present at <40 [iM in the upper lay- 
ers of Garden Banks sediments. Their concentrations were, how- 
ever, below our detection limits of 0.5 u,M in the sediment slurries 
at both To and Tf. Even if one assumes that 40 u,M NO^ was 
present in sediments, it could only produce 50 u,M DIC (5 u,mol 
C) through complete NO ^ reduction to N2. Thus, nitrate 
coupled alkane oxidation cannot solely explain the observed 
discrepancies. 

To determine if alkanes were systemically lost via leakage 
and/or other bottle effects, sediment-free reactors were incubated 
in parallel, and alkane concentrations were monitored over the 
course of the incubations. These reactors exhibited <5% loss over 
the course of the incubations, which is markedly lower than our 
least active biological treatment (CH4). SRR is known to be sensi- 
tive to the incubation time, although sediments — including those 
from the GoM — are typically incubated for 24 h as was used here 
(Fossing and Jorgensen, 1989). Moreover, the rates shown here 
are comparable to previously published rates (Arvidson et al., 
2004; Joye et al., 2004). In any case, shorter incubation times 
might result in higher SRR, and might account for some of the 
observed excess alkane consumption. Leakage of gas from sam- 
ple vials (different than the incubation vials) could also affect the 
isotopic signatures of the alkane pools and could be a source of 
error resulting in lower precision measurements (though nothing 
in our data is consistent with this hypothesis). 

Shifts in the microbial community could lead to enrichment of 
acetoclastic methanogens that can use acetate (a possible product 
of partial propane or butane degradation) to produce methane. 
Methane has been reported as a potential carbon sink dur- 
ing degradation of higher molecular weight hydrocarbons (Gray 
et al., 2010). However, in the present study we do not observe any 
changes in the concentration of the methane produced compared 
to the N2 controls (Table S2). We also do not observe an increase 
in the total methanogen population in the C2-C4 amended sedi- 
ments by qPCR with primers specific for the methanogenic mcrA 
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gene. Thus, methane appears unlikely to be a carbon sink in our 
experiments. 

Finally, precipitation of carbonates may be a sink of carbon 
within the sediments. The inorganic precipitation of carbonate, 
known as authigenic carbonate, can occur at the sediment-water 
interface or within the sediment pore water. Authigenic carbon- 
ates are often formed in sediments where increasing alkalinity, 
typically from sulfate or metal reduction, increases the carbonate 
saturation state past a saturation threshold causing precipitation 
of minerals (calcite or aragonite). Sulfate-reduction increases the 
alkalinity of pore waters by removing hydrogen ion from the 
local environment in the form of H2S and generating bicarbonate 
concentrations by oxidizing organic carbon. Authigenic carbon- 
ates are found throughout the GoM (Roberts and Aharon, 1994; 
Sassen et al., 2004). Generally, environments such as the GoM 
with substantial amounts of organic carbon but that hinder aero- 
bic respiration and support alkalinity-increasing processes such as 
sulfate reduction have the potential to harbor large carbon sinks 
in the form of authigenic carbonate (Higgins et al., 2009). While 
inorganic precipitation of carbonates is possible in our incuba- 
tions, data for the total sedimentary inorganic carbon content is 
unavailable. 

Microbial community analyses via pyrotag sequencing impli- 
cate that certain members of the class Deltaproteobacteria are 
enriched during the batch incubations. Further phylogenetic 
analyses indicate that the enriched bacteria are closely related to 
previously enriched/isolated C3-C4 degrading SRP (Kniemeyer 
et al, 2007; Jaekel et al., 2013), as well as uncultured marine 
SRP observed in GoM sediments (Orcutt et al, 2010; Kleindienst 
et al., 2012). The Deltaproteobacterial sequences most enriched 
in the ethane incubations were closely related to isolate BuS5, 
and the enrichment culture Butanel2-GMe, both of which belong 
to the DSS cluster (Kniemeyer et al, 2007). Intriguingly, previ- 
ous studies suggest that BuS5 (degrades propane and n-butane) 
and enrichment Butanel2-GMe (degrades n-butane) do not 
degrade ethane (Kniemeyer et al, 2007; Jaekel et al., 2013). 
Other uncultured DSS cluster members were also identified in 
the ethane degrading incubations (Kleindienst et al., 2012). Thus, 
it is possible that SRP closely related to the C3-C4 degrad- 
ing DSS cluster might be associated with ethane degradation 
in these incubations, though this hypothesis remains to be 
tested. 

While pyrotag sequencing using 16S rRNA gene data show 
the phylogenetic structure of the microbial community (and is 
not quantitative), qPCR analysis helps us to quantitatively assess 
the functional potential of the microbes in the reactors. Many of 
the Deltaproteobaceria reduce sulfate using the dissimilatory sul- 
fite reductase or adenosine 5'-phosphosulfate reductase enzymes 
(encoded by dsrA and aprA genes, respectively). The qPCR results 
(both dsrA and aprA gene abundance) demonstrate that the total 
number of SRP decrease compared to the To sample and the 
N2 control. However, C2-C4 consumption correlates stoichio- 
metrically with SR, and SR rates were higher in the presence of 
the alkane gases compared to the N2 control. Collectively, these 
results indicate that there maybe an enrichment of a specific sub- 
set of the SRP community responsible for the consumption of 
C2-C4. While lower in abundance than the total SRP in the Tq 



sample or N2 control, this subset of the community likely exhibits 
higher specific SR activity. 

The data presented here provide insight into alkane oxida- 
tion rates in Garden Banks sediments. Like previous studies, these 
data confirm that the addition of alkanes stimulates sulfate reduc- 
tion (Figure 2, Table 4). Notably, the rates of C2-C4 consumption 
are comparable to CH4 consumption, though their stoichiomet- 
ric impacts on the sulfate pool vary. For example, assuming 
complete oxidation, it is likely that 1, 1.75, 2.5, and 2.8 moles 
of sulfate are reduced per mole alkane for C1-C4, respectively. 
Consequently, the relative effect of C2-C4 oxidation on the sul- 
fate pool is much greater than for methane given the observed 
similarity in the oxidation rates. Accordingly, alkane oxidation 
may represent a substantial sink for sulfate in sediments where 
alkanes are elevated, such as the GoM where they can constitute 
more than 10% of the total gas pool (Milkov, 2005). This is most 
likely relevant deeper in the sediments, where alkane concentra- 
tions are highest and where sulfate concentrations are lowest. In 
such scenarios, it is not implausible that C2-C4 oxidation might 
limit the availability of sulfate for methane oxidation, though 
this speculation requires further study. The data further constrain 
carbon exchange between the alkane and DIC pool, and this phe- 
nomenon should be considered when interpreting DIC isotope 
ratios of alkane-replete sediments. 

Ethane is the next most abundant short-chain non-methane 
alkane at our study site (Table 3). However, previous stud- 
ies have reported ethane-driven sulfate reduction at very slow 
rates by microbial enrichments obtained from a similar loca- 
tion (Kniemeyer et al., 2007). Anaerobic ethane oxidation likely 
involves a novel mechanism because it requires the activa- 
tion of a primary carbon, in contrast to butane where sec- 
ondary carbons are available (Kniemeyer et al, 2007, and 
references therein). Notably, we observed substantial ethane 
consumption over the course of these incubations [approxi- 
mately two orders of magnitude higher than those reported 
by Kniemeyer et al. (2007), and comparable to the oxida- 
tion rates of methane, propane or butane]. The presence of 
ethane also stimulated sulfate reduction, which implies a rela- 
tionship between these processes. Our data show anaerobic 
ethane utilization begins without delay, similar to other alka- 
nes, suggesting that microorganisms are poised for ethane 
oxidation. 

Bulk geochemical and isotopic surveys of alkanes along seeps 
have been used to imply microbial consumption of short-chain 
alkanes (Sassen et al., 2004; Orcutt et al., 2005). Our data confirm 
that short-chain alkanes are oxidized to DIC, likely coupled to sul- 
fate reduction. However, both our alkane consumption rates and 
the isotopic shifts observed in the DIC pools suggest other sinks 
may exist. At thermogenic hydrocarbon seeps, these processes 
may have an important impact on the local carbon and sulfur 
cycles. The strategies used herein — namely the combination of 
molecular, geochemical, and isotopic assessments — was leveraged 
to establish the relationships between anaerobic alkane oxidation 
and SRRs, carbon flux, microbial activity, and microbial com- 
munity composition and phylogeny. Future experiments should 
consider these and previous data to gain further insight into the 
signatures and mechanisms of these biogeochemical processes 
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as well as organisms involved in anaerobic short-chain alkane 
oxidation. 
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